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Abstract 

The bound state energies and wave functions for a particle exposed to the Hulthen potential field 
in the D-dimensional space are obtained within the improved quantization rule for any arbitrary 
I state. The present approximation scheme used to deal with the centrifugal term in the effective 
Hulthen potential is systematic and accurate. The solutions for the three-dimensional (D = 3) 
case and the s-wave (/ = 0) case are briefly discussed. 
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I. INTRODUCTION 



It is well-known that the exact analytical solution of the hyperradial Schrodinger equation, 
in any arbitrary spatial dimension (D > 2), for its bound state energy levels, is fundamental 
in understanding the bound energy spectrum of nonrelativistic and relativistic quantum me- 
chanics, since the resulting wave function contains all the necessary information to describe 
fully any quantum system. There are only a few potentials for which the radial Schrodinger 
equation can be solved explicitly for all n and I quantum numbers. One of these exactly 
solvable potentials is the Hulthen potential [1,2] which can be solved in a closed form for s 
wave (1 = 0). However, the three-dimensional radial Schrodinger equation for the spherically 
symmetric Hulthen potential cannot be solved analytically for I ^ states because of the 
centrifugal term ~ r~ 2 [3-5]. The Hulthen potential is one of the important molecular poten- 
tials used in different areas of physics to describe the interaction between two atoms and has 
attracted a great of interest for some decades in the history of quantum chemistry. Until now, 
it has also been used extensively to describe the molecular structure and a possible form of 
atomic and nuclear interactions [6-10]. Further, the study of this potential is only a special 
example of those exponential type potentials [11] such as the modified hyperbolic- type po- 
tentials (Scarf, modified Rosen- Morse, modified second type Poschl- Teller) [12], Eckart [13], 
Rosen- Morse [14], Manning- Rosen [15] potentials and so forth. So far, numerous attempts 
have been developed to calculate the bound-state energies such as variational [10], supersym- 
metry quantum mechanics [16,17], Nikiforov-Uvarov method (NU) [5,14,18,19], asymptotic 
iteration method (AIM) [20], hypervirial perturbation [21], shifted 1/N expansion (SE) and 
modified shifted 1/N expansion (MSE) [22], exact (improved) quantization rule (EQR or 
IQR) [23,24], perturbative formalism [25-29], polynomial solution [30], wave function ansatz 
method [31], factorization method [32] and tridiagonal J-matrix representation (TJM) [33] 
which split original Hamiltonian into two parts as H = H + V where H Q is the part of the 
Hamiltonian that could be treated analytically while the remaining part, V, has to be treated 
numericallyto solve the radial Schrodinger, Klein-Gordon and Dirac equations exactly or 
quasi-exactly for I ^ within a given potential. 

Recently, Ma and Xu have proposed an exact (improved) quantization rule (EQR or 
IQR) and shown its power in calculating the energy levels of all bound states for some 
solvable quantum systems [23,24]. The method has been shown to be effective for calculating 
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the bound state solutions of the Schrodinger and Dirac wave equations with a spherically 
symmetric potential [34-39]. So far, it has been applied, with great success, to study a great 
number of potentials like the rotating Morse [34,35], the Kratzer-type [36], the trigonometric 
Rosen-Morse [37], the hyperbolic and the second Poschl- Teller-like potentialas [38], the 
Hulthen potential [39] and the Woods-Saxon potential [40] and so forth. Very recently, Gu 
and Sun [39] have extended the application of the IQR to the solution of the .D- dimensional 
Schrodinger equation with the Hulthen potential for I ^ using the usual approximation to 
deal with the centrifugal term [41-43]. 

Very recently, Dong [12] has introduced a more beautiful exact quantization rule to 
simplify the calculation of the energy levels for exactly solvable quantum systems. The 
energy spectra of the modified hyperbolic-type potentials have been carried out by this 
rule. Qiang and Dong [44] have found a proper quantization rule (PQR) and showed that 
the previous complicated and tedious calculations for the energy spectra can be greatly 
simplified. This new quantization rule can be applied to any exactly solvable potential. 
Qiang-Dong PQR has been further applied to exactly solvable shape invariant potentials 
[45]. Very recently, Yin et al. [46] have shown that the SWKB is exact for all shape 
invariant potentials (SIPs). 

In this paper, we aim to extend the study of Ref. [39] by using a new improved ap- 
proximate scheme to deal with the centrifugal term. Further, we solve the present potential 
model on the assumption that the space may posses an arbitrary number of spatial dimen- 
sions D. This arbitrary dimensional study enables one to give analytical tests using energy 
calculations for interdimensional degeneracy, i.e., (n, l,D) — > (n,l±l,D=p2) corresponding 
to the confined D = 2 — 4 dimensional Hulthen potential. 

It is worth noting that this alternative approximating approach has shown its accuracy 
in calculating the analytic and numerical energy spectrum of the Hulthen potential for / ^ 
[3-5]. Further, it has also been applied to the spin and pseudospin symmetries, e.g., Wei and 
Dong have studied the approximation of the Dirac equation with scalar and vector modified 
and deformed generalized Poschl- Teller and Manning-Rosen potentials within the improved 
approximation formula to the centrifugal term [47-50]. 

This paper is organized as follows. In Sec. 2, the EQR (IQR) method is reviewed and 
extended to any arbitrary dimension (D > 2). In Sec. 3, the D-dimensional (D > 2) 
Schrodinger equation is solved by this method with I ^ quantum numbers to obtain the 
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energy eigenvalues. In Sec. 4, we calculate the corresponding hyperradial wave functions of 
the Hulthen potential. Finally, some conclusions are given in Sec. 5. 

II. EXACT (IMPROVED) QUANTIZATION RULE 

A brief outline to the improved quantization rule is presented with an extension to the 
D-dimensional space (D > 2). The details can be found in Refs. [23,24]. The IQR has 
recently been proposed to solve exactly the one-dimensional (ID) Schrodinger equation: 

r(x) + k(xf^(x) = o, k( X ) = vgg^E ZM, (i) 

where the prime denotes the derivative with respect to the variable x. Here \i is the reduced 
mass of the two interacting particles, k(x) is the momentum and V(x) is a piecewise con- 
tinuous real potential function of x. The Schrodinger equation is equivalent to the Riccati 
equation 

<t>'(x) + (/)(x) 2 + k(xf = 0, (2) 

where <p(x) = ip' (x) / ip (x) is the logarithmic derivative of wave function ip(x). Due to the 
Sturm-Liouville theorem, the <f>(x) decreases monotonically with respect to x between two 
turning points, where E > V(x). Specifically, as x increases across a node of the wave 
function ip(x), <p(x) decreases to — oo, jumps to +oo, and then decreases again. 

Moreover, Ma and Xu [23,24] have generalized this exact quantization rule to the three- 
dimensional (3D) radial Schrodinger equation with spherically symmetric potential by sim- 
ply making the replacements x — > r and V(x) — > V c s(r): 

fB r B 

fk{r)dr = N« + fk' { r)f±dr, k(r) = vgggZE^W l, (3 ) 

where r A and r B are two turning points determined from the relation E n j = V e s(r), N = n+1 
is the number of nodes of <p(r) in the region E n ^\ > V c g(r) and it is larger by one than the 
number of nodes of wave function ip(r). The first term Ntt is the contribution from the 
nodes of the logarithmic derivative of wave function, and the second term in (3) is called 
the quantum correction. It is found that, for all well-known exactly solvable quantum 
systems, this quantum correction is independent of the number of nodes of wave function 
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of the system. This means that it is enough to consider the ground state in calculating the 
quantum correction, i.e., 

rB 

f k \r)^p-dr = n q , (4) 
J 0o ( r ) 

The quantization rule still holds for Schrodinger equation with spherically symmetric po- 
tential in D dimensions. In what follows, we shall employ this method to extend the work 
of Ref. [39] by using an improved approximation to the centrifugal term. 

III. EIGENVALUES OF THE HULTHEN POTENTIAL 

The Schrodinger equation with spherically symmetric potential V(r) for I ^ takes the 
simple form 

(~h v2 ° + v{r) ~ En '*) ^'' m(r ' n ° ] = °' (5) 

where the representation of the Laplacian operator V|,, in spherical coordinates, is 



D dr 2 r dr r 2 

and 



x? 



1>nj, m (r, tl D ) = ^,(0 = r-< D -WR(r), (7) 

where Y^C^d) is the hyperspherical harmonics. The wave functions ip nlm (r,Qi)) belong to 
the energy eigenvalues E n i and V(r) stands for the Hulthen potential in the configuration 

/ D \ V2 

space and r represents the D-dimensional intermolecular distance j . 

\i=i 

Further, substituting Eqs. (6) and (7) into Eq. (5) yields the wave equation satisfying 
the radial wave function R(r) in a simple analogy to the 2D and 3D radial Schrodinger 
equation 

R "( r ) + p [ E n,i ~ V e ff(r)\ R(r) = 0, (8) 

where V e ff(r) is the Hulthen effective potential in D dimensions defined by 

. . o e~ ar (A 2 - 1) h 2 
V ef Ar) = -Ze 2 a- + v - ! — , 9 

with the parameter 

A = 2/ + D-2. (10) 
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The radial wave function R(r) satisfying Eq. (8) should be normalizable and finite near r = 
and r — > oo for the bound-state solutions. The wave equation (8) with the Hulthen potential 
is an exactly solvable problem for I = (s-wave) [2,51-53], however, it cannot be solved 
analytically when I ^ because of the centrifugal barrier term, i.e., (A 2 — 1) h 2 r~ 2 /(8fx). 
Therefore, to solve Eq. (8) analytically, we must use a new approximation scheme of the 
exponential- type proposed recently by Jia et al (cf. e.g., [54-59]) to deal with the centrifugal 
term: 

1 / e~ ar \ 

3~« 2 c o + - — 2 , (11) 



r z \ (1-e 

where the dimensionless constant cq = 1/12 is exact as reported by other authors (cf. e.g., [3- 
5]). Very recently, we have applied the above approximation scheme (11) to obtain improved 
bound state solutions to the Schrodinger equation with the Manning-Rosen potential for 
arbitrary /-waves [15]. Obviously, the above approximation to the centrifugal term turns to 
r~ 2 when the parameter a goes to zero (small screening parameter a) as 



lim 

a->0 



2 1 * 1 
a I c + — 7 + ^— — 



' (Ha) 



..2 



which shows that the usual approximation is the limit of our approximation (cf. e.g., [4] 
and the references therein). Further, by defining 

a = Ze 2 a, b = c?L\ L 2 = ^(l + ^) (l + ^) , (12) 



then we have from Eq. (8): 

2/i 



R"(r) + 



h 2 



End + a- — - b c + 



R{r) = 0, (13) 



1 - e~ ar \ (1 - e - c 
where E n j is the bound state energy of the system and n and I signify the radial and angular 
quantum numbers, respectively. 

We now study this system through the improved exact quantization rule. At first, we 
introduce a new variable 

e -ctr 

<r) = 1 _ e . ar , z'(r) = -az(l + z), (14) 

where r G (0, oo) and z G (0, oo). Overmore, the turning points z\ and zb are determined 
by solving V e ff(z) = bz 2 + (b — a) z + bc = E n j as follows: 



z a =Yb~l~^b^ a ~ b)2 + 4b ( Enl ~ bco) ' ZB = Tb~l + Tb^ a ~ b)2 + 4b (E ™ " 6co) ' 

(15) 



with the properties 



<2 E n l tirx\ 

za + z B = ^ - 1, z A z B = — + cq. (16) 



The momentum fc(z) between two turning points is expressed as 

/ (a -b) (bc - £„,/) R 

k{z) = — — y -z 2 + — - — z = —^—y/(z B -z)(z-z A ), (17a) 



ft 



<i/c(z) \/2/i6 / [zb — z I z — za 



dz 2h \ \ z — za V z b — z 
The Riccati equation (2) now becomes 



;i?b) 



_ az {z + i) = _^ [ Eo _ ^ + (Q _ b) z _ 6 i _ (jz)2j (18) 

dz n 

having the only possible solution satisfying 

o (r) = ciz + c 2 , o '( r ) = -olcxz{\ + z), ci > 0. (19) 

where we have used o ( r ) = ^ol- 2 )- Substituting (p (z) into Eq. (18), one has the ground 
state energy eigenvalue and wave function solutions 



o (z)=^ + a(^-f), 



"=5 + W 1 + ^ = ' + ^ 1 . ( 2 °) 



2 1 2y 1 r 1 2 ' — ' 

h n=0 -ho-bco- — {-J^^ ~ 2 ) 



After a lengthy algebra but straightforward, we can calculate the integral of the quantum 
correction (4) based on the ground state as 



Q c = nq = n(^L + v-lj. (21) 
The integral of the momentum k(r) in the quantization rule (3) is calculated as 

'''B ZB / , , 

J f,^ dr = _ J ( V( Z ~ Z A ) {ZB - Z) _ y/(z - Z A ) (Zb ~ z) j ^ 



Using the relations (21) and (22), the improved quantization rule (3) turn out to be 



k—L 1 + yco-- r -yc + - h = 7T I —A + n + u). (23) 



Thus, one can finally find the approximation to the bound state energy levels E n i for the 
D-dimensional Hulthen potential, 



E 



, D-3' 
I + I c - 



fiZe 2 



h 2 {n + l + V) 



a 



(24) 



where n, I = 0, 1, 2, • • • . Therefore, the energy spectrum in 3D space can be obtained as 



E. 



n,l 



2/z I 12 



fiZe 2 



[n + l + 1) 



-i 2 



ft 2 (n + Z + 1) a 



, n,Z = 0,l,2, 



(25) 



which is identical to Eq. (34) of Ref. [4]. In the case of the s-wave (/ = 0), the previous 
relation turns out to become 



En 



h a 
'~2fT 



fiZe 2 



(n + l) 



-i 2 



h 2 (n + 1) a 



, n = 0,l,2, 



(26) 



which is identical to the ones obtained before using the AIM [20], SUSYQM approach 
[60-64], quasi-linearization method [65] and NU method [42,66]. Overmore, if we take the 
dimensionless constant c = in the present approximation, Eq. (24) reduces to 

2 



2^,2 



E n l 



h a 
~2/7 



jiZe 



(n + l + ^1) 



h 2 (n + l + a 



(27) 



which is consistent with the energy eigenvalues formula given in Eq. (32) of Ref. [20], Eq. 
(24) of Ref. [67] and Eq. (28) of Ref. [42] when D = 3. By taking the chosen parameters 
h — 2/j, — e — 1 and for Z — 1, the above result is consistent with Eq. (24) of Ref. [39]. 
The critical screening parameter can be found as a c = h 2^^ +1 ^ when E n \ = and c = 0. 

IV. EIGENFUNCTIONS 



We are now in the position to study the corresponding eigenfunction of this quantum 
system for completeness. The Riccati equation of the relation (8) is 



where 



Hr) = -^[E nl -V eff (r)]-ci>(ry, 



R'(r) 



R(r) 



(28) 



(29) 
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Based on 

R(r) = e* r * {r)dr = e~ - f (30) 
and using Eq. (19), we can easily calculate the eigenfunction of the ground state as 

R (r) = N (e- ar f (l - e~ ar )\ e > 0, v > 1, (31) 

where 

/ 2/1 iiZe 2 1 v . . 

V ft cr n va 2 

with v is defined in Eq. (20) and N is the normalization constant. 

Let us find the eigenfunction for any quantum number n. At first, considering the bound- 
ary conditions 

{0 when r — >■ oo, 
(33) 
1 when r — > 0, 

with — >■ 0, based on Eq. (31), we may define a more general radial eigenfunctions, valid 
for any quantum number n, of the form: 

R(y)=y^{l-yYF(y), y = e~ ar , e n , t > 0, !/ > 1, (34) 

satisfying the boundary conditions in Eq. (33), where 

fiZe 2 1 n + u 

£ n,l = . 2 ( ~ v > 0. 35 

Substituting Eq. (34) into Eq. (8) leads to the following hypergeometric equation 



y(l-y) F"{y) + [1 + 2? n ,, - (1 + 2e n)l + 2v) y] F'(y) - 
whose solutions are the hypergeometric functions 



+ 2e nt i) - 



2fiZe 2 



h 2 a 



F{y) = 0, 
(36) 



k=0 

where 



jp( \ v (A nr \ T(C) ^ T(A + k)T(B + k) y k 
F(y) = 2Fl (A, B; C; y) = ^ f ^ T ^ W (37) 



,_ 2 2fiZe 2 
A = e n j + u- \/e nt i + = -n, 



B = e n , l + v+ x lel l + '^^, (38) 



C = 1 + 2e 



n.l- 
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By considering the finiteness of the solutions, the quantum condition is given by 



e n ,i + u- d£ 2 n j + j^- = - n , n = 0, 1, 2, • • • , (39) 

from which we obtain Eq. (25). Now, we may write down the radial wave functions (34) as 

R(r) = N nl (e- ar Y n > 1 (1 - e~ ar Y 2 F, (-n, n + 2 (e^ + + 2e nX , e~ ar ) . (40) 

If we set n = in Eq. (40), then we can easily obtain Eq. (31). Finally, the unnormalized 
total wave functions are obtained as 

i> n ,U r ' Q d) = Kir- (D - 1)/2 (e- ar Y n ' 1 (1 - e~ ar Y 2Fl (-n, n+2 (e B>I + v) ; l+2e B>I ; e~ ar )YT(n D ). 

(41) 

which is identical to Eq. (42) of Ref. [4] when D = 3. Thus, the Jacobi polynomials can be 
expressed in terms of the hypergeometric functions [69] 

Pi A > B \l - 2s) = r ^|/ + + ^ 2F1 n + A + B + 1; A + 1; *). (42) 

The hypergeometric function 2-P\(A, -B; C;s) is a special case of the generalized hypergeo- 
metric function [69,70] 

Man « 2 , • • • , op; • • • , /3 S ; x) = jr ^Hf^T^S' ( 43 ) 

where the Pochhammer symbol is defined by (y) k = T(y + k)/T(y). 

Let us find the normalization constant. Introducing the change of parameters y{r) = 
e~ ar and making use of Eq. (41), with the help of Eq. (42), we are able to express the 
normalization condition J °° i?(r) 2 (ir = 1 as 



a 



[ y 2 ~^-\l - y) 21 ^- 1 [P^ 2l+D - 2 \l - 2y)] 2 dy = 1. (44) 
Jo 



Unfortunately, there is no formula available to calculate this key integration. Neveretheless, 
we can find the explicit normalization constant N n[ . For this purpose, it is not difficult to 
obtain the results of the above integral by using the following formulas [69-71], 

n 1 / r — 1 \ n ~ p f r 4- 1 \ p 

P^( X ) = { n + a)l { n + fll g ^ - - _ ?y (/? - ?y ^ - ^ (— ) (_), 

(45) 
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and 



B 



(x,y)= f ^{l-tf-'dt 
Jo 



V(x)F(y) 



T(x + y) 

Thus, the normalization constant M n i is now obtained as 



, Re(x),Re(y) > 0. 



Mr 



nl 



(n + 2l + D- 2)\T(2e n j + n + 1) 



i 



ar(2e ra ,i + 27i + 2f + D + l) 

n 

T(2e n>l + 2n + l) £ (f P f q f M ) 

p,q=0 



-1 



where 



(46) 



(47) 



f p = {-iy P \T{2Z n + + n - p + 1)(2Z + p + D - 2)\ (n + p)\, (48a) 

/, = (-l) q q\r(2e nl + n-q+l)(2l + D + q-2)\(n + q)\, (48b) 

f p , q = (2l+p + q + D-l)L (48c) 

It is worth noting that one of the disadvantages of the EQR (IQR) approach is that it 
cannot get the eigenfunctions of studied potential models. The estimation shown in this 
section is only the reverse of logarithmic derivative to the original Schrodinger equation, 



i.e., (p (r) 



dr 



ln(i? ( r )) [ see Eq. (30)] [72]. The traditional method is thus used again. 



V. CONCLUSIONS 



In this work, we have applied an alternative method to obtain approximate energy eigen- 
values and eigenfunctions of the D-dimensional Schrodinger equation for the Hulthen po- 
tential with I 7^ within the improved approximation scheme for the centrifugal term. The 
advantage of this method is that it gives the eigenvalues through the calculation of two 
integrations (21) and (22) and solving the resulting algebraic equation. First, we can eas- 
ily obtain the quantum correction by only considering the solution of the ground state of 
quantum system since it is independent of the number of nodes of wave function for exactly 
solvable quantum system. Second, the wave functions have also been obtained by solving 
the Riccati equation. The general expressions obtained for the energy eigenvalues and wave 
functions can be easily reduced to the three-dimensional space (D — 3), s-wave (/ = 0), the 
Co = (usual approximation) cases. The method presented here is a systematic one, simple, 
practical. and powerful than the other known methods. It is worth to extend this method 
to the solutions of other nonrelativistic [23,24,34,36-39] and relativistic [35] wave equations 
with different potential fields. Finally, it can be also used to deal with many exactly solvable 
quantum systems with wide range of potentials as stated by many authors [23,24,34-39]. 
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Appendix A: Integral Formulas 

The following integral formulas are useful during the calculation of the momentum integral 
and the quantum correction terms [34,36,68]: 

rB 




(49) 




(50) 




(51) 




(52) 




(a + br) yj(r - r A )(r B - r) 



1 



dr 



y/(a + br A )(a + br B ) 



7T 



r B > r A > 0. 



(53) 
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